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A simple, but fully-covariant model for describing neutral Vector Mesons, in both light and heavy sectors, is 
briefly illustrated. The main ingredients of our relativistic constituent model are i) an Ansatz for the Bethe- 
Salpeter vertex for Vector Mesons, and ii) a Mandelstam-like formula for the electromagnetic decay widths. 
The free parameters of our approach are flxed through a comparison with the valence transverse-momentum 
distribution, n{kj_), obtained within phenomenological Light-Front Hamiltonian Dynamics models reproducing 
the mass spectra. Preliminary results for the transverse-momentum distributions, the parton distribution and the 
electromagnetic decay constants are shown. 



1. INTRODUCTION 

Aim of this contribution is to illustrate a fully 
covariant model for describing the electromag- 
netic (em) decay of neutral Vector Mesons (VM) , 
in both light and heavy sectors (see [1] for a 
preliminary presentation of our approach). For 
the present stage, we have considered only the 
ground states, to be described by using Bethe- 
Salpeter (BS) amplitudes with a simple analytic 
form, namely exhibiting only single poles. This 
assumption allows us to easily perform the ana- 
lytic integration needed for calculating the con- 
stituent quark (CQ) transverse-momentum dis- 
tribution inside the VMs (see, e.g. Refs. [2] and 
[3] for the pion case) and then em decay widths. 
The free parameters present in the Ansatz for 
each VM are fixed by comparing the so-called va- 
lence transverse-momentum distribution, n{k±), 
with the same quantity evaluated within a Light- 
Front Hamiltonian Dynamics (LFHD) approach 
exploiting the eigenfunctions of two phenomeno- 
logical mass operators, presented in Refs. f4| and 



O, and able to reproduce the VM spectra. It 
should be pointed out that the free parameters 
represent an effective way for including some non 
perturbative features in our analytical Ansatz. 

A possible form of the BS amplitude for an in- 
teracting qq system with J = 1, can be written 
as follows 

*A(fc, P) S{k, mi) Vx{P) S{k ~ P, m^) x 

AyAf(fc,F) (1) 

where 5(fc, m) — S{kon,m) +7+/(2fc+) is the 
Dirac propagator of a constituent quark with 
mass m, fc^„ = {fcQ„,fc+,kj_} with fc~^ — {m? -\- 
|k±P)/A;+, P'^ the four-momentum of a VM with 
mass = Myj^,j, Kvuik^P) the momentum 
dependence of the BS amplitude, and Vx{P) = 
£\{P) ■ V{P) with e^(P) the polarization four- 
vector (A is the helicity) and V'^{k,k — P) the 
Dirac structure given by the following familiar ex- 
pression, transverse to P'^, (see, e.g., 6J 
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Then, the scalar product Vx{P) reduces to 

h[MvM-f] 
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It should be pointed out that one could have two 
different AyM(fe, -P)'s, one for each Dirac struc- 
ture in V\{P) (i.e. /a and ), but at the present 
stage we assume a simpler form, with the same 
AvM{k,P) for the two structures, that leads to 
the expected Melosh Rotation factor for a sys- 
tem, in the limit of non interacting systems, i.e. 
MvM Mfree Or P'^ P^/^^^ 0. Therefore the 
following comparison with the phenomenological 
LFHD outcomes can be made more strict. 

In the preliminary calculations presented in 
this contribution, the momentum dependence 
of the BS amplitude is given by the following 
Ansatz, with only single poles, viz 
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[fc2 -ral + {P~ kf 
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where m^;. , i = 1, 2, 3 are free parameters (deter- 
mined as described below), and M the normaliza- 
tion factor, that is obtained by imposing the nor- 
malization for the BS amplitude in Impulse Ap- 
proximation, i.e. by adopting constituent quark 
free propagators. The form chosen for Kynik, P) 
allows one: i) to implement the correct symme- 
try under the exchange of the quark momenta 
(for equal mass constituents), as follows from the 
charge conjugation of the neutral mesons; ii) to 
regularize the integrals needed in our approach 
for the evaluation of valence wave functions, 
decay constants, transverse-momentum distribu- 
tions, etc.; and iii) to avoid any free propagation 
in the valence wave function, given the presence 
of the numerator in AvA/(fc, P) (see also |ll8j ). 

To determine m^j. in Eq. ((4]), we have intro- 
duced the transverse-momentum distribution of a 
CQ inside the VM, adopting the same definition 
already applied by de Melo et al. [5] for a co- 
variant description of the pion. In a frame where 
= 0, one first defines the valence component 



of the VM state in the standard way (see, e.g., 
0), viz. 

^'^t,(e,kx;mHj = j dk- S{kon,mi) x 

Vx{P) AvAiik, P)S {{k - P)on, ma) (5) 

with ^ — k'^ / P^ . The normalization factor 
in KvM{k,P) (cf Eq. @) is evaluated in Im- 
pulse Approximation by exploiting the total- 
momentum sum rule, that reads for the plus com- 
ponent as follows 

[ AvM ik,P) 
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In Eq. JT]), T+{P,k,on) and T+{P,k,ist) are 
the on-shell and the instantaneous contributions, 
respectively, according to the well-known de- 
composition of the Dirac propagator (see be- 
low Eq. (HJ). Their explicit expressions will 
be given elsewhere [TU]. Notice that, after per- 
forming the trace, the dependence upon the he- 
licity disappears, and therefore its presence in 
TJ_s^(P, fc) becomes dummy. Finally, let us in- 
troduce the valence transverse-momentum distri- 
bution, n„o/(fc_L), given by 
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ponent, given by 
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where P^^ is the probability of the valence com- 
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Table 1 

The mass eigenvalues for the ground states of the neutral VMs investigated in this contribution compared 
with the corresponding experimental values. The mass operator corresponds to the one in Ref. [3] (see 
[TU] for details) . The quark masses adopted for the actual evaluation are also listed. 



VM 


mq{MeV) 


M'y\j (MeV) 


M'yZ (MeV) 


P 


328 


777.0 


775.500 ± 0.4 




449 


1020.7 


1019.455 ± 0.020 




1559 


3082.9 


3096.916 ± 0.011 


T 


4891 


9537.2 


9460.300 ± 0.268 



Table 2 

The same as in Table 1, but for the mass operator of Ref. [5]. 



VM 


mq{MeV) 


M'y\, (MeV) 


M'yZ (MeV) 


P 


220 


777 


775.500 ± 0.4 




419 


1016 


1019.455 ± 0.020 




1628 


3091 


3096.916 ± 0.011 


T 


4977 


9460 


9460.300 ± 0.268 



rfki T+(P,k,on) |$^tf(^,ki;™flJI' (9) 



and 



$rM(e,k^;™i?J= fc+ (p-fc)+ j dk- 

AvM{k,P) 

(P-m?) [(fc_p)2_^2] 



(10) 



In a LFHD approach (see, e.g. [2] and refer- 
ences quoted therein), one has 



|^yM(e,kx)P (11) 

where iPvm{£,i^i) is the eigenfunction of a rel- 
ativistic mass operator and Mj^^^ — (m^ + 
|k_Lp)/^(l — ^). In order to implement our pro- 
cedure for fixing the free parameters in Avm, 
we have first evaluated ny^j{k±) by adopting 
the eigenfunctions obtained by using the model 
mass operators of both Refs. and jj], and 
then we have minimized the difference | Uyai {k±) — 
ny^{k±)\. It should be pointed out that the VM 
spectra are satisfactorily reproduced for each VM 
investigated in this contribution. As an example 



of the reliability of the adopted LFHD models, 
the masses of the considered neutral VMs for the 
ground states are compared with the correspond- 
ing experimental values in Tables [T] and [5] for the 
mass operators of Refs. [4j and [5j, respectively. 
Indeed, also the masses of the first three/four (de- 
pending upon the considered VM) excited states 
are well reproduced by the model mass operators. 

Once the parameters have been determined, 
one can calculate both static (e.g., em decay 
rates) and dynamical (e.g. parton distributions) 
quantities. For completing this section, devoted 
to the general presentation of the formalism, we 
give the expression of the parton distribution, i.e. 



(12) 



qiO = / dk^ fci /i(C,fcj.) 



in terms of the chiral-even transverse-momentum 
distribution, fi{^,k±), (see, e.g. [10] and [3] 
for the pion case), that yields the distribution 
of a constituent inside the meson with the lon- 
gitudinal and transverse components of the LF- 
momentum. It reads 



/i(^,fci) = 
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The explicit expression of fi{C,k±) will be pre- 
sented elsewhere [T^. 

2. THE MANDELSTAM FORMULA 
FOR THE EM DECAY CONSTANT 

In order to evaluate the em decay constant, 
fvM, we adopt a Mandelstam-like formula 
(see also [H] and [3] for the pion case). The start- 
ing point is the macroscopic definition of fvMi 
through the transition matrix clement of the em 
current for a given neutral VM, viz 



(0|J^(0)|F,A) =zV2/vMe^ 



(14) 



Let us remind that the decay constant fvM is 
related to the em decay width as follows 



^ e+e- 
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In our model, the transition matrix element in 
Eq. (|14p can be approximated microscopically a 
la Mandelstam through 

AvM{k,P,mi,m2) ^ 

Tr[VA(P) + m2)r{^ + "^i)] (16) 
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where 



2p = 



(Qu — Qd) 

n/2 



S,//* — Qc Qt — Qb 

with Qi the quark charge. The lengthy algebra for 
evaluating the analytic integrations in Eq. (jl6|) 
will be reported elsewhere [10]. 

3. PRELIMINARY RESULTS 

In Figs. 1 and 2, the valence transverse- 
momentum distributions (cf Eq. ([5])) for p, </>, 
J/ip and T are shown, for different sets of the 



free parameters. In order to get some dynamical 
input, in Fig. 1 the free parameters of our Ansatz 
(see Eq. ([1])) have been determined through the 
model of Ref. while in Fig. 2 through the 
model of Ref. [5 . In the figures the distributions 
are multiplied by a factor fcj^ , to have an immedi- 
ate intuition of the transverse-momentum region 
relevant in the evaluation of the valence probabil- 
ities. It is worth noting that the long tail of the 
heavy VM is an expected feature. In Figs. 3 and 
4, the parton distributions are presented. As it is 
shown, the general behavior at the end-points and 
the accumulation of the distribution for a; = 1/2 
as the CQ mass increases are recovered by using 
our Ansatz. 

Finally, in Tables 3 and 4, the preliminary val- 
ues for the em decay widths, for the two different 
choices of the free parameters, are compared with 
the corresponding experimental data. It should 
be pointed out that while the model of Ref. [3] 
yields reasonable results for the light sector (let 
us remind that the confining interaction for the 
squared mass operator of ^ is parabolic), in the 
heavy sector the model of Ref. seems to behave 
better (for this model a linear confining interac- 
tion is considered in the mass operator). Such 
a comparison suggests that an improved descrip- 
tion of the confining interaction for the squared 
mass operator of Ref. [4] could lead the theo- 
retical predictions of the em decays closer to the 
experimental ones. 

4. CONCLUSIONS 

The preliminary results for some static and 
dynamical quantities for neutral VM, obtained 
within our covariant description, both in the light 
and the heavy sectors, have been shortly pre- 
sented. After calculating the valence wave func- 
tions, Eq. ([5]), and the transverse- momentum 
distributions, Eq. ([8]), we have fixed the 
three parameters in our Ansatz for the Bethe- 
Salpeter amplitude through a comparison with 
the transverse-momentum distribution obtained 
within a Light-Front Hamiltonian Dynamics ap- 
proach. The eigenvectors of two different mass 
operators, corresponding to Ref. [4 and Ref. [5], 
have been considered in order to perform the com- 
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Table 3 

Preliminary VM em decay widths, calculated using our Ansatz, ([1]), and the free parameters determined 
through the comparison with the eigenfunctions of the model 4 for the mass operator. 



VM 


rf+,- (keV) 


r:::- (kev) m 


P 


8.579 


7.04 ± 0.06 




1.952 


1.27 ± 0.04 


Jli, 


2.526 


5.55 ± 0.14 


T 


0.187 


1.340± 0.018 



Table 4 

The same as in table 1, but using the eigenfunctions of the model [5] to fix the free parameters. 



VM 


r*'|,- (keV) 


VZl- (keV) m 


P 


24.791 


7.04 ± 0.06 




4.342 


1.27 ± 0.04 


Jli, 


7.102 


5.55 ± 0.14 


T 


0.493 


1.340± 0.018 
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Figure 1. Valence transverse-momentum distri- 
butions for a constituent inside p, J/^* and 
T vs the quark transverse momentum. The nor- 
malization is J dk± nyai{k±) = 1. Solid line: 
p; dashed line: </>; dotted line: J/^*; double-dot- 
dashed line: T. The free parameters in ([T]), have 
been determined through the comparison with 
the eigenfunctions of the model [4j . 



Figure 2. The same as in Fig. 1, but using the 
eigenfunctions of the model [5j to fix the free pa- 
rameters in ([T]). 
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0.8 1 



parison. Then, we have evaluated the parton dis- 
tributions and the em decay constants (cf Table 
Eland Table SI). 

The work in progress will substantially improve 
the present calculations, in two respect: both in- 
troducing a more refined Ansatz for the BS ampli- 
tude and considering new phenonienological mass 
operators. 
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